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Abstract 



We compare different models for hadronic and quark phases of cold baryon-rich matter 
in an attempt to find a deconfinement phase transition between them. For the hadronic 
phase we consider Walecka-type mean-field models which describe well the nuclear sat- 
^f) ' uration properties. We also use the variational chain model which takes into account 

correlation effects. For the quark phase we consider the MIT bag model, the Nambu- 
, Jona-Lasinio and the massive quasiparticle models. By comparing pressure as a func- 

tion of baryon chemical potential we find that crossings of hadronic and quark branches 
are possible only in some exceptional cases while for most realistic parameter sets these 
^ ■ branches do not cross at all. Moreover, the chiral phase transition, often discussed within 

the framework of QCD motivated models, lies in the region where the quark phases are 
unstable with respect to the hadronic phase. We discuss possible physical consequences 
of these findings. 

X 

1 Introduction 

It is commonly believed that quarks and gluons are relevant degrees of freedom in strongly 
interacting matter at very high temperatures and baryon densities. This state of matter is 
usually described by various QCD motivated models. On the other hand, at low temperatures 
and moderate densities, at least up to the nuclear saturation density p Q = 0.17 fm~ 3 , strongly 
interacting matter exists in the hadronic phase. In particular, atomic nuclei are finite droplets of 
this phase with baryon density ps — Po- They are self-bound and therefore can exist in vacuum, 
without external pressure. This fact itself provides an important constraint on the equation 
of state (EOS) of cold nuclear matter, namely, its pressure must vanish at ps = Po- There 
are many effective models which successfully describe nuclear matter in terms of interacting 
nucleons. 



Unfortunately, at present there exists no rigorous approach which can describe the EOS 
of strongly interacting matter at finite baryon densities. As known, QCD lattice simulations 
have principal limitations at nonzero chemical potential. Therefore the only practical way to 
study the possibility of the deconfinement phase transition in this case is to compare various 
models of hadronic and quark phases. In the last two decades there were numerous attempts 
to construct a unified EOS which would interpolate between the two asymptotic regimes. In 
this paper we critically revise this problem in the light of new calculations for the hadronic and 
quark phases. 

In our analysis we include four hadronic models. Two of them are Relativistic Mean-Field 
(RMF) models of the Walecka type |1|, §, namely, the NLZ and the TM1 § models. Next is 
the so called Chiral Hadronic Model (CHM) which was recently developed in Refs. []5], |6|. Fourth 
model is the Variational Chain Model (VCM) based on the Argonne NN potentials with 
addition of 3-body forces and relativistic corrections. Unlike the above three models, the VCM 
takes into account correlation effects neglected in the mean-field approximation. As argued in 
Ref. J7|, the mean-field approximation is not well justified at baryon densities p# < p • 

For the quark matter we take three different models: the MIT Bag Model (BM) ||, the 
Massive Quasiparticle Model (MQM) @, [10| and the Nambu-Jona-Lasinio model (NJL) 0, [TJ]. 
Presumably, the hadronic degrees of freedom are more relevant at baryon densities ps < Po 
and quark degrees of freedom take over at much larger densities. But the transition between 
these two regimes is very poorly understood at present. By considering various hadronic and 
quark models we pursue several goals. First, we compare different models for a single phase to 
get an idea on the uncertainty in their predictions. Second, by applying the Gibbs criterion, 
we investigate the possibility of a deconfinement phase transition in cold baryon-rich matter. 
Finally, we examine reliability of different quark models (in particular, allowed values of model 
parameters) by extrapolating their predictions into the domain of nuclear matter, p# ~ p . 

The paper is organized as follows. In Sect. II we give short descriptions of several popular 
models of the hadronic phase. Their predictions regarding the equation of state of cold baryonic 



matter are summarized and compared between each other. In Sect. Ill different models of the 
quark phase are introduced. Their pairwise comparison with hadronic models, aimed at finding 
a hadron-quark phase transition, is carried out in a systematic way. In Sect. IV we present 
our conclusions and outlook. 

2 Models of hadronic phase 
2.1 Relativist ic mean— field models 

At present the field-theoretical description of dense hadronic matter is one of the most popular 
approaches. Within this approach the matter is described in terms of baryons interacting with 
self-consistent meson fields. Most calculations are done within the mean-field approximation. 
There are many versions of the RMF model which differ by the choice of meson fields as well as 
by the baryon-meson coupling schemes. Here we consider two realizations of the RMF approach 
which give very good description of finite nuclei, namely, the NLZ || and the TM1 |3| models. 
The general form of the effective RMF Lagrangian used in these models is[] 

C = ifj (i$ - m* N - s^7o) i> - U s {a) + U v (u) . (1) 

Here if), a and u are, respectively, the nucleon, scalar and vector meson fields, 

m* N = m N - g a a (2) 

is the effective nucleon mass, U s and U v are the scalar and vector potentials: 

Us(a) = Mr^ + ^V + ^V* (3) 

U v (u) = + . (4) 

In the above equations m$ denote vacuum masses of nucleons (i = N) and mesons (i = a,u), 
gj are coupling constants. The parameter sets for the NLZ and TM1 models are listed in 

Table Note that contrary to the NLZ, within the TM1 model the vector field is a nonlinear 

1 Below we consider static and homogeneous isospin-symmetric matter. In this case the derivatives of mean 
meson fields over space and time as well as the contribution of p-mesons may be omitted. 



function of the baryon density. This leads to a significant reduction of the repulsive interaction 
at high p B . 

Table 1: Parameters of RFM models 





ttin (GeV) 


m CT (GeV) 


m^(GeV) 


9a 


92 (fm" 1 ) 


93 


9u 


93u 


NLZ § 
TMl H 


938.9 
938.0 


488.67 
511.198 


780 
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10.0553 
10.0289 


-13.5072 
-7.2325 


-40.2243 
0.6183 


12.9086 
12.6139 




71.305 



Within the mean-field approximation a and u fields are regarded as purely classical and 
replaced by c-numbers. Using the Lagrangian ([]]) one can easily calculate pressure P of cold 
nonstrange matter as a function of the baryon chemical potential ps (see for details Refs. |], 0). 
At given \ib , applying thermodynamic relations one obtains the following equations for the 
baryon density and the energy density of matter 

p B = dP/dp B , (5) 

e = PbPb - P • (6) 

The energy per baryon is equal to E/B = e/ps ■ 

2.2 The chiral hadron model 

Although the Lagrangian (|l|) leads to very good description of nuclear phenomenology it has 
one principal defect. Namely, it does not respect chiral symmetry of strong interaction. It is 
commonly accepted that this symmetry is spontaneously broken in vacuum and can be restored 
at high density and temperature. In recent years there were several attempts to incorporate this 
symmetry into the RMF framework. It turned out that most easily this can be achieved [|T3"| , p~4"f 
by introducing an additional scalar (dilaton) field x responsible for the trace anomaly of QCD. 
For our analysis we have chosen the Chiral Hadron Model (CHM) developed in Ref. ]5|, [J. 
Below we use the version of the CHM, given by the parameter set CI in Ref. ||. As compared 
to the NLZ and TMl models, the CHM includes also a strange scalar field (. The following 
parametrization of m* N is used instead of (@): 

™*n = 9a° + 9c ■ (7) 



The scalar potential U s is parametrized as a fourth-order polynomial in a, ( and x- Its param- 
eters are tuned to reproduce the pattern of spontaneous symmetry breaking at low densities. 
In addition there appear logarithmic terms, oc x 4 l°gX 4 an d oc x 4 l°g(cx 2 £), motivated by the 



and finite nuclei. 

2.3 The variational chain model 

The VCM is a microscopic approach where the interaction between nucleons is describes in 
terms of two and three body forces. The Hamiltonian is written in the form 



Here the first term gives non-relativistic kinetic energy of nucleons. The NN interaction 
terms are chosen in the form of the Argonne potential with parameters fitted to reproduce 
the NN scattering phases up to the cm. energy 300 MeV. The terms Svij denote relativistic 
(Lorentz-boost) corrections up to the second order in the NN pair momentum. The 3N inter- 
action terms Vijk are taken in the Urbana UIX form with parameters which correctly reproduce 
binding energies of lightest (A < 4) nuclei. The calculations are carried out by applying the 
variational Monte Carlo method and the chain summation technique [^]. Besides properties of 
finite nuclei, pressure and energy density of nuclear matter are calculated at various nucleon 
densities and neutron to proton ratios. 

2.4 Comparison of model predictions 

In Figs. we compare pressure and binding energy of cold isospin-symmetric baryonic matter 
calculated within the above four models. Figure [l] shows predictions for energy per baryon as 
a function of ps ■ All models give nearly the same results at densities ps < Po • At higher 
densities the CHM, VCM, and the TM1 model are not far from each other, unlike the NLZ 

model which seems to overestimate significantly the energy per baryon at ps > 3poQ- 

2 One should regard the results of effective hadronic models at ps » Po with caution due to the lack of 
available information concerning EOS at such densities. 



trace anomaly. As shown in Ref. || the CHM gives satisfactory description of nuclear matter 
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Figure 1: Energy per baryon of cold nuclear matter calculated within different hadronic models 
as function of baryon density. Here and below the curves are specified in the key box. 




Figure 2: Pressure of cold nuclear matter vs baryon chemical potential calculated within dif- 
ferent hadronic models. 



Figure |2| shows pressure as a function of baryon chemical potential. The discrepancy between 
different hadronic models becomes especially evident at pB—m N > 200 MeV. One can see that 
the results of the NLZ model strongly deviate from the predictions of other models. Apparently, 
this is due to overestimation of vector repulsion within the NLZ model. On the other hand, the 
CHM, VCM, and the TM1 model give similar results at \i& < 1.5 GeV. All hadronic models 
predict the first order phase transition of the liquid-gas type at ps — itln which corresponds 
to baryon densities ps < Po- This is clear from Fig. || where pressure is shown at a finer scale. 
For example, in the case of the VCM the parts AB, BC, and AC of the pressure curve are 
unstable with respect to decomposition of matter into the dilute {i = 1) and dense (i = 2) 
phases. According to the Gibbs rules for coexisting phases 

p(i) = P (2) j (9 ) 

/ig> = a4 2) - (10) 

The phase transition points at T = are given by intersection of different branches of pressure 
as a function of chemical potential 0. For example, point A in Fig. ^| is the phase transition 
point in the case of the VCM. 

According to Eq. (H), the density jump in the first order phase transition, — Pb is equal 
to the difference of slopes dP/dps at pb — * p c ± where p c (the phase transition point) is the 
solution of Eq. (|9]). The equivalent criterion of the first order phase transition is the so-called 
"double tangent construction" for energy density as a function of baryon density. Indeed, by 
using thermodynamic identities (§)-(@) one can rewrite Eqs. ©-([H]) in the form 

eg) - e W _ deU _ de^ 
P ( b ~ Pb dpB dpB 
Below the conditions ([|)-(|ll]) and ( [TT]) will be used for investigating the possibility and esti- 



mating parameters of the deconfinement transition by comparing pressures of hadronic (i = 1) 

and quark (i = 2) phases. 

3 When analyzing EOS in the /is — P plane one should bear in mind that lower branches of pressure 
correspond to metastable or unstable states. 




Figure 3: The same as in Fig. ^ but for small values of baryon chemical potential. Points A,B, 
and C show boundaries of metastable and unstable regions in the case of the VCM. 
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Figure 4: Baryon density of cold nuclear matter as function of chemical potential. 
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Figure 5: Pressure vs baryon density calculated within different hadronic models. 

Four hadronic models considered in this paper predict essentially different behavior of the 
baryon density with raising ps ■ This is shown in Fig. |j. Again, one can see strong deviation 
of the NLZ results from the predictions of other hadronic models. Figure |5| shows pressure 
isotherms as functions of p#. It is clearly seen that the NLZ model predicts the "hardest" 
EOS. Indeed, at intermediate densities ps ~ (1.5 — 4) p this model predicts the highest com- 
pressibility i.e. the largest derivative dP/dps- The NLZ model gives also the largest pressure 
at given ps ■ On the other hand, the P{pb) curve predicted by this model is lower as compared 
to other models (see Fig. ^|). Apparently, this difference is caused by the vector repulsion which 
is especially large in the NLZ model. 

An important characteristic of nuclear matter is the in-medium nucleon mass m* N . Ex- 
periments on inelastic electron-nucleus scattering show that at ps ~ po , this mass is lowered 
to about 0.7rri]y . By solving the gap equations @, (0) one can calculate the p# dependence 
of m* N within the CHM as well as in the TM1 and NLZ models. The results are shown in 
Fig. |]. All the models predict a rapid drop of m* N at p# < 2p and much slower decrease at 
higher densities. This behavior shows the tendency towards restoration of chiral symmetry at 
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Figure 6: Nucleoli effective mass as function of baryon density predicted by different hadronic 



models. The dots show estimates [HI obtained from the QCD sum rules. 



high p B . 

Having in mind that predictions of the NLZ model at high densities strongly deviate from 
the results of other models, below we use only the CHM, VCM, and the TM1 model for the 
comparison with quark models. 

3 Deconfined phase 
3.1 The bag model 

The simplest description of deconfined phase is given by the MIT bag model (BM). Within this 
model pressure of homogeneous nonstrange quark matter at zero temperature is expressed as 



P = P (fi)-B, 

where \i = /is/3 is the chemical potential of light quarks, B is the bag constant and 
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Figure 7: Comparison of energy per baryon of cold symmetric matter calculated within the 
MIT bag and different hadronic models. Values of S 1//4 in MeV are given in parentheses. 

is pressure of ideal gas of massless quarks with Nf flavors. Below we assume that only light 
u, d quarks are present, i.e. Nf = 2. We consider the simplest version of the BM neglecting 
perturbative corrections due to the color-magnetic interaction. Thus, the only nonperturbative 
aspect of the model is associated with the bag constant B. The BM has been widely used for 
modelling a phase transition from the baryon-free hadron matter into the quark-gluon plasma 
at nonzero temperature. The critical temperature T c for this transition is determined not only 
by the bag constant but also by the numbers of active degrees of freedom in respective phases. 



Reasonable values, T c ~ 160 MeV, are obtained [16| with B 1 ^ ~ 200 MeV. The calculations 



below are made for J B 1 / 4 =145, 165 and 200 MeV. 

In Figs. |7]-|9]we compare the EOS predicted by the hadronic models (VCM, CHM and TM1) 
and the BM for different values of B . Figure [7] shows energy per baryon as a function of baryon 
density. It is clear that small values of bag constant, _B 1//4 < 150 MeV, are unrealistic because 
in this case quark matter would be more stable at ps ~ Po than cold nuclear matter. Due to 
this reason, below we consider only two values of bag constant corresponding to B 1 ^ = 165 




M B (GeV) 

Figure 8: The same as in Fig. [7j but for pressure vs baryon chemical potential. Dots (triangles) 
mark the points of possible phase transitions from the VCM (TM1) hadronic phase to the BM 
quark phase with B 1 / 4 = 165 and 200 MeV. 

and 200 MeV. 

Figure |8| shows pressure as a function of chemical potential calculated for the same models. 
From this figure one can see that at B l l A = 165 and 200 MeV the transition from the hadronic 
to the quark phase is possible only for the VCM and the TM1 model. Filled symbols in 
Fig. |8| mark points of intersection between respective pressure curves. By using the Gibbs 
rules (|9|)-(jl0"D it is possible to calculate parameters of the quark-hadron mixed phase predicted 



Table 2: Parameters of hadron-quark phase transition predicted by matching the hadronic 
models and the BM with different bag constants B . 
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by these models. The results of such calculation are summarized in Table |2|. A more detailed 
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Figure 9: The same as in Fig. [7|, but for pressure vs baryon density. Pairs of dots (triangles) 
connected by dashed lines show mixed hadron-quark states obtained by matching the EOS, 
predicted by the VCM (TM1 model) and the BM with B l ' A = 165 and 200 MeV. 

information on the resulting EOS with the hadron-quark phase transition is given by Fig. || 
However, the calculations show that the phase transition occurs at rather large values of fiB 
and pb, in the region where applicability of hadronic models is questionable. 

3.2 The massive quasiparticle model 



A more realistic model of the quark phase, the MQM, has been proposed in Refs. ||, [T0| 



The case of hot baryon-free matter has been studied in Ref. |J. In Ref. [|Hj a similar model 

was formulated to study the EOS of quark matter at T = 0, \i ^ 00. Within the MQM it 

is assumed that cold quark matter may be regarded as an ideal gas of quarks with nonzero 

effective mass m = m(/i). Arguments in favor of this picture follow from calculations based 

on the hard thermal loop resummation technique developed in Refs. |T7], [18], In the case 

of hot baryon-free matter it was possible to reproduce lattice results by using only leading 

order diagrams for quark and gluon self-energies. On the other hand, at small strong coupling 
4 It should be noted that the dependence of strong coupling constant on \i was disregarded in Ref. Jul . 



constant, a s = g 2 /An , one obtains results consistent with the perturbation theory up to the a 2 
order. The case of cold quark matter has been recently studied |5(J within the hard density loop 
resummation technique. It was shown that dependence of pressure on quark chemical potential 
can be well approximated by treating quarks as ideal gas of massive quasiparticles. However, 



the analytic formula P = P^{m^) suggested in Ref. [^0| is thermodynamically inconsistent in 
the case when m(/z) ^ const (see below). 

As shown in Refs. [[H], ^2], |23f , at large enough 3-momenta the dispersion relation for quarks 
can be interpreted in terms of quasiparticles with nonzero effective mass m . In the leading 
order in ot s this mass can be expressed as |^3| 



m 



2a, 
3n 



(14) 



One should have in mind that at small values of /i corresponding to a s > 1 (see below) the 
applicability of Eq. (|i~4D becomes questionable. 



In our calculations we use the three-loop expression [24] for the strong coupling constant 
a s = g 2 /An as function of the renormalization scale Q 
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Here L = \og(Q 2 /A 2 ) and 



A, = 11 - 2JV//3, Pi = 51 - 19AT//3, 
p 2 = 2857 - 5033^/9 + 325JV//27 . 



(16) 



The cutoff momentum A is fixed by the condition a s (2 GeV) = 0.3089 [P4| , which gives 
A ~ 0.4178 GeV for Nf = 2. From dimensionality arguments it is clear that Q ~ fi . As 
in Refs. |I| ^0|, p5]| , it is assumed that Q = 7/i where the coefficient 7 is of the order of unity 
and does not depend on fi . At fixed 7 physical values of /i correspond to Q > A . 

Further it is postulated that the density of quark matter equals to the density of ideal gas 
of massive fermions 
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Figure 10: The QCD coupling constant as function of baryon chemical potential at Nf = 2 
and different values of the parameter 7 . Thin dashed line shows maximal possible value of a s 
within the MQM. 

where Q(x) = (1 + sgnx)/2. Thus, the quark density vanishes at fi < m. This condition 
is satisfied at \i < fi c where \x c is the critical chemical potential determined from Eq. (|b|) 
with m = fi . The latter condition is equivalent to a s = 3n/2. The calculation shows that 
/i c ~ 0.5974/7 GeV for N f = 2. 

Figs. [10|-|n| show a s and m as functions of baryon chemical potential /ie = 3/z for several 
values of the parameter 7. The shaded region in Fig. |TT] corresponds to the domain where 
m > fi . It is seen that quark masses are very sensitive to the choice of 7. At large hb quark 
masses increase, although slower than /!# . As discussed in Ref. [f26j, the existence of nonzero 
quark mass does not contradict to the restoration of chiral symmetry at large /is . 

Using Eq. (|T?| ) and thermodynamic identity (|5]) one can calculate pressure as a function 
of \i : 



Here we assume that -P(/i c 



P = J dfixpaifn) . (18) 

He 

, i.e. pressure is zero at vanishing quark density. If one 
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Figure 11: Effective mass of u, d quarks m within the MQM vs baryon potential of quark 
matter \i B . Unphysical region is shown by shading. 

would replace m(/ii) in the integrand pid by the "constant" mass m(p) , one would obtain the 
well-known expression for pressure of ideal gas Pid- However, as shown in Ref. [27|, the 



approximation P = P^ is thermodynamically inconsistent when masses depend on T or fiB ■ 
In particular, Eq. (|5|) does not hold in this case. The thermodynamic consistency can be 
recovered by introducing the effective bag constant P(p) = P;a — P . 

Figure [12] shows thermodynamic quantities calculated within the MQM by using Eqs. fll7|)- 
(|T8|), (P). Density, pressure and energy density are shown as functions of the dimension- 
less variable /i//i c . All quantities are given as ratios to their respective limits for massless 
quarks (Psb = £sb/3 = Po(fi)). In this representation the results do not depend on 7 . In 
Fig. [T2| we also demonstrate inaccuracy of the approximation P = used in Ref. PU| . It is 
seen that P(/x) is typically of order O.lpd and changes sign at \i ~ 2/i c . 

Energy per baryon and pressure calculated within the MQM are compared with predictions 



of hadronic models in Figs. |i3Hl4 Unlike the BM, in this model the energy per baryon increases 




Figure 12: Pressure P, energy density e and number density p of quarks within the MQM (nor- 
malized to their Stephan-Boltzmann limits, m — > 0) as functions of quark chemical potential \i 
in units of its critical value fi c . Thin solid and dashed lines show pressure and energy density 
calculated by using the ideal gas approximation B(p) = 0. 
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Figure 13: Comparison of energy per baryon of cold symmetric matter calculated within the 
MQM and different hadronic models. 



1000 



800 - 



600 - 



VCM 

TM1 

CHM 

— mqm(t=i) 

MQM (7=1.5) 

MQM ( 7 =2) 



> 
CD 




400 - 



200 - 



Figure 14: The same as in Fig. [13|, but for pressure vs baryon chemical potential. Filled symbols 
show points of possible phase transitions. 

monotonically with p#, even at low densities. According to Fig. [13], the MQM with 7 > 2 is 
unrealistic since in this case the quark matter would be energetically more favorable at ps % Po 
than normal nuclear matter. 

According to Fig. |H| at 1.5 < 7 < 2 the MQM pressure curve intersects only with curves 



predicted by the VCM and the TM1 model. Parameters of possible phase transitions are given 
in Table [3[ One can see that quark phase appears only at rather high baryonic densities, namely, 
at pb > Pb > 6po- At 7 < 1 there are no intersections of hadronic and quark branches, at 
least at not too large ps- 

Table 3: Parameters of quark-hadron phase transition predicted by matching the VCM, 
TM1 model, and MQM at different 7 . 
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We have also checked sensitivity of the results to the choice of Nf in the QCD coupling 
constant. The calculation has been made using Nf = 3 in Eqs. (|T3j)-(r^), but retaining Nf = 2 
(only u and d quarks) in thermodynamic quantities. This leads to a slight increase of P(/j>b) 
that lowers the phase transition density by not more than 10%. 

3.3 The Nambu— Jona-Lasinio model 



The NJL model JL1], [12| is one of the most popular models dealing with constituent quarks. 
There are several advantages of the NJL model compared to other models considered above. 
First, it respects chiral symmetry of strong interactions, second, it explicitly takes into account 
negative energy (Dirac sea) states, and third, it describes well meson phenomenology. Different 
versions of this model have been extensively used to describe the EOS of equilibrium and 
nonequilibrium quark matter at finite T and fiB ■ 

Here we use results of our calculations [fl^, |29|, [3(J within the SU(3)-flavor version of the 



model suggested in Ref. ||31|| , but with an additional term due to vector-axial- vector interaction. 



The color-singlet part of the Lagrangian in the mean field approximation can be written as 

f f 

+ ^Uph + ^WPsf (19) 
/ / 

Here ipf is the field operator of quarks with flavor f = u,d,s and 

ps f = <^ f ip f >, (20) 
p V f = < ^floipf > (21) 

are their scalar and vector densities. Angular brackets in Eqs. (PD|)-(PTD denote quantum- 
statistical averaging. Gs, Gy and K in Eq. fliTf ) are, respectively, the coupling constants of 
scalar, vector and flavor-mixing interactions. 

The constituent quark masses, mf, are determined from the coupled set of gap equations 

rrif = m of - G s Psf + 2K Yl ps f ' ( 22 ) 



where mo/ is the bare (current) mass of quarks with flavor / . 

The NJL model is an effective, non-renormalizable model. To regularize the divergent 
contribution of negative energy states of the Dirac sea, one must introduce an ultraviolet 
cut-off. Below the 3-momentum cut-off G(A — p) is used in divergent integrals. The model 
parameters m f,Gs,K,A can be fixed by reproducing the observed masses of tt,K , and 7/ 



mesons as well as the pion decay constant. As shown in Ref. plj , a reasonable fit is achieved 
with the following input parameters: 

m 0u = m od = 5.5 MeV, m 0s = 140.7 MeV, (23) 
G s = 20.23 GeV~ 2 , A = 0.6023 GeV, K = 155.9 GeV" 5 . (24) 

In principle, the vector coupling constant may be extracted by fitting the nucleon axial 
charge or masses of vector mesons. It was shown that the ratio £ = Gy/Gs should be of 



the order of unity. However, as discussed in Ref. P3" |, the accuracy of such fitting procedure 
is rather low. Due to uncertainty in the parameter Gy , below we present results for various 
values of £ from the interval < £ < 1 . 

Let us consider cold isospin-symmetric nonstrange matter with baryon density ps • In this 
case py s = 0,py u = pyd = 3ps/2. The energy density e can be calculated directly from the 
Lagrangian ([19]) with the result [p9] 



E eKf + ^Psf+^fPvf ~4^n^/ + e °- ( 25 ) 



/ 



Here ejcf is the kinetic term which includes also "active" negative energy states with mo- 
menta p < A: 



PFf 

3 



e Kf = — I dpp x /mf+p 2 , (26) 



A 



where PFf = ( 71 " 2 Py//3) 1,/3 is the /-quark Fermi-momentum. Scalar densities can be calculated 
by using the relation 

\dm f ) 



PFf 
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Figure 15: Comparison of energy per baryon of cold symmetric matter calculated within the 
MQM and NJL model. 



The constant eo in the l.h.s. of Eq. fl25|) is introduced in order to set the energy density of the 
physical vacuum equal to zero. This constant can be calculated by solving the gap equations (P^) 
for the vacuum case, ps = . Finally, pressure can be found by using thermodynamic identities 
Pb = de/dp B and P = PbPb — e • 

The results of calculations are shown in Figs. [E)|-|T%| for three values (0, 0.5, 1.0) of the 



parameter £ . Figs. |l5HL6l show the comparison of the MQM and NJL model. It is seen that 
the discrepancy between their results is especially large in the case 7 = 1. Unlike the MQM, 
the NJL model predicts metastable states of quark matter corresponding to local minima in 
E/B as function of ps . However, such states exist only at £ < 0.10. 



In Figs. 17-18 we compare EOS predicted by the NJL and hadronic models. As seen in 



Fig. [18] the hadron-quark phase transition may take place only at small £ which, most likely, are 
not realistic. For larger £ the repulsive vector interaction makes the NJL phase too stiff to cross 
any of the hadronic curves. The phase transition parameters for cj = are given in Table HI 



As shown in Rcfs. 1291, 30 , they appear in strange matter even at larger £. 
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Figure 16: The same as in Fig. [15], but for pressure vs baryon chemical potential. 



Again the transition to quark matter is possible only from the TM1 and VCM hadronic phases 
and predicted critical densities are above 5po- 

At £ < 0.71 the pressure isotherms calculated within the NJL model contain unstable parts 
in the region \±b — 1 — 1.2 GeV |29j . This means that this model itself predicts the first-order 
chiral phase transition in the quark matter. However, corresponding parts of pressure isotherms 
lie below the hadronic curves. Therefore, such a phase transition is not observable due to the 
hadronization of quark phase. This situation is quite general: many phase transitions found 
within different quark models are predicted in regions of the \ib ~ P plane where quark phase 
is unstable with respect to hadronization. 



Table 4: Parameters of hadron-quark phase transition matching the VCM and TM1 model 
with the NJL model (f = 0). 
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Figure 17: Comparison of energy per baryon of cold symmetric matter calculated within the 
NJL model and different hadronic models. 
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Figure 18: The same as in Fig. [D], but for pressure vs baryon chemical potential. Filled symbols 
show points of possible phase transitions. 



4 Conclusions and outlook 



Our goal in this paper was to find possible phase transition between hadronic and quark phases 
of strongly interacting matter at high baryon density. For each of the phases we have selected 
several models which are widely discussed in the literature. Namely, we have considered the 
CHM, VCM, and RFM models for the hadronic phase, and the BM, MQM, and NJL model 
for the quark phase. When necessary, we have varied model parameters in reasonable limits 
to determine possible range of uncertainty. For each model we have calculated main thermo- 
dynamic quantities, in particular, pressure as a function of baryon chemical potential, P(hb) ■ 
Possibility of a deconfinement phase transition is signalled by an intersection of P{hb) curves 
corresponding to the hadronic and quark phases. 

A pairwise comparison of different hadronic and quark models has shown that such an 
intersection is not at all a general feature but rather an exception. For instance, the CHM, 
which is probably the most advanced hadronic model at present, does not predict a phase 
transition with any of quark models. Also, no phase transition is predicted when most realistic 
parameters are used in quark models (B 1 ^ ~ 200 MeV in the BM, 7 ~ 1 in the MQM, or 
£ > 0.5 in the NJL model). Even in the cases when a phase transition is possible, e.g. between 
the TM1 model and the BM, or between the VCM and BM, characteristics of mixed phase are 
very sensitive to model parameters. Moreover, this phase transition is predicted at such a high 
density, above 5po where predictions of hadronic models are very unreliable. The situation 
with quark models is even worse. First, it is unclear at all how far down in density one can 
use these models. Second, the predictions of various quark models differ significantly in the 
region of moderately high baryon densities of interest here. Of course, due to the asymptotic 
freedom of QCD the quark phase must approach asymptotically the ideal gas limit. However, 
it is unclear at present when such asymptotic behavior sets in. 

In this paper we have considered only isospin-symmetric matter, but we think that similar 
conclusions can be made also for /3-equilibrium neutron star matter. For example, the com- 
parison of EOS predicted by the CHM and NJL model (for details see Ref. |33|) does not show 



any phase transition between hadronic and quark phases. 

As has been demonstrated in this paper, some phase transitions found in the QCD motivated 
models, e.g. the chiral transition within the NJL model, occupy regions of the //# — P plane 
where quark phases are unstable with respect to hadronization. We expect that similar situation 



takes place also for color-superconducting phases of quark matter |34|, |35|, although this question 
deserves a special study. 

We should conclude that after more than 30 years of model building in both the hadronic 
and quark sectors the situation regarding a hadron-quark phase transition and EOS at high 
baryon densities remains rather uncertain. We believe that progress in this field can be achieved 
by developing new class of models where both the hadronic and quark degrees of freedom are 
treated within a unified theoretical framework. Attempts to construct such a model were made 



in Refs. |36| , T7[| . Such unified approach should include constraints from nuclear physics (exis- 
tence of the nuclear bound state) and QCD (chiral symmetry, asymptotic freedom). It is quite 
possible that the transition from hadronic to quark degrees of freedom will be continuous 



like ionization in atomic systems [38[ . Indications of such behavior are found in recent lattice 



calculations | 39| . All this means that the problem of hadron-quark transition at high baryon 



densities remains a challenge for theorists. 
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